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Abstract 

The fragmentation of quasi-projectiles from the nuclear reaction ^°Ca+^^C 
at 25 MeV/nucleon was used to produce excited states candidates to a- 
particle condensation. Complete kinematic characterization of individual 
decay events, made possible by a high-granularity 47r charged particle multi- 
detector, reveals that 7.5±4.0 % of the particle decays of the Hoyle state 
correspond to direct decays in three equal-energy a-particles. 
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Bose-Einstein condensation is known to occur in weakly and strongly in- 
teracting systems such as dilute atomic gases and liquid "^He yy]. During the 
last decade it was theoretically shown that dilute symmetric nuclear mat- 
ter may also experience Bose particle condensation jl, Is], More precisely, 
for densities smaller than one fifth of the nuclear saturation density, nuclear 
matter organizes itself in a-clusters, while for higher densities the 2-nucleon 
deuteron condensation is preferred. This new possible phase of nuclear mat- 
ter may have its counterpart in low-density states of self conjugate lighter 
nuclei, in the same way as superfiuid nuclei are the finite-size counterpart 
of superfiuid nuclear and neutron matter. This means that under some cir- 
cumstances, the alpha condensation, i.e. bosonic properties, might dominate 
over the nucleon properties even in finite nuclei. 

Considerable advance in favour of a-particle condensation in nuclei is pro- 
vided by the excellent theoretical description of the Hoyle state {i.e. the first 
0+ excited state at 7.654 MeV of ^^C) and of the state at 15.097 MeV of 
^^O in terms of condensate type wave functions [s], [oj. The case of the Hoyle 
state is particularly suggestive as both shell-model and no-core shell model 
calculations are known to fail in describing it. Recent Fermionic Molecular 
Dynamics (FMD) calculations j3] plead in favour of a more nuanced interpre- 
tation of the Hoyle state structure, that is a mixture of various pre-formed 
a-configurations. This scenario agrees with the pioneering works of Uegaki 
et al. [8| where acute-angle, bent and linear chain configurations were iden- 
tified. In what regards the expected diluteness, coupled channel analysis 
of experimental data confirm theoretical works jsj estimating the rms radius 
of the Hoyle state as 45% larger than the radius of ^^C in its ground state. 
More generally, a-particle condensation is conjectured to be a generic feature 
of medium-size self-conjugated 4iV nuclei whose excitation lies in the vicinity 



of the Na decay threshold [10|, lU 



The aim of the present Letter is to search from the experimental side a 
direct evidence for a-particle condensation from the Hoyle state. Accord- 
ing to the present understanding of the Hoyle state: a gas-like structure of 
three a-particles which occupy dominantly the lowest S orbit, such an ev- 
idence may be judged from the simultaneous emission of three a-particles 
with very low kinetic energy dispersion. As it is experimentally difficult to 
identify and measure a-particles with low energy, probably the most appro- 
priate strategy, chosen in this work, should involve high velocity reaction 
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products in the laboratory (to take advantage of velocity boosts) detected 
by a high granularity-high solid angle particle array (to precisely reconstruct 
the directions of velocity vectors). The existence of the direct decay channel 



was firstly investigated in Ref. [12| and its contribution to the total alpha 
decay (dominated by the ^Bec,.^. + a sequential decay) was estimated to be 
at most 4%. The existence of such a decay channel is also of the greatest 
importance for nuclear astrophysics, due to the crucial role played by the 
Hoyle state for the synthesis of ^^C in the universe 13|, |l4 . 

The data reported here have been obtained in the nuclear reaction '^'^Ca+^^C 
at 25 MeV per nucleon incident energy performed at INFN, Laboratori 
Nazionali del Sud in Catania, Italy. The beam impinging on a thin carbon 
target (320 /xg/cm^) was delivered by the Superconducting Cyclotron and the 
charged reaction products were detected by the CHIMERA 47r multi-detector 



151]. The beam intensity was kept around 10'' ions/s to avoid pile- up events. 
CHIMERA consists of 1192 silicon-CsI(Tl) telescopes mounted on 35 rings 
covering 94% of the solid angle, with polar angle ranging from 1° to 176°. 
The solid angle corresponding to each module is not uniform, but varies be- 
tween 0.13 msr at forward angles and 35.4 msr at the most backward angles. 
Among the most interesting characteristics of CHIMERA are the low detec- 
tion and identification thresholds for light charged particles (LCP) and the 
very high granularity at forward angles. The mass and charge of the detected 
nuclei were determined by the energy-time of flight method (TOP) for LCP 
stopped in silicon detectors and AE — E {Z > 5) and shape identification 
{Z < 5) techniques for charged products stopped in CsI(Tl). In addition 
^Be nuclei (two equal-energy a's hitting the same crystal) were identified in 
CsI(Tl) [16]. The energy of detected nuclei was measured by the Si detectors 
calibrated using proton, carbon and oxygen beams at various energies rang- 
ing from 10 to 100 MeV. Por Z = 2, dedicated energy calibrations of the fast 
component of CsI(Tl) light was realized using the TOP. Though high-quality 
TOP charts allowing for a direct calibration of the Csl light exist only for 
60% of the total number of modules, we have finally managed to calibrate 
more than 95% of modules from 1° to 62°. The modules for which TOP 
information was poor or missing were calibrated by comparing the fast com- 
ponent distribution with the benchmark distribution of the corresponding 
ring built out of the telescopes with excellent TOP. The energy resolution 
for alpha particles varies between 1.0 and 2.5% depending on the module. 

Complete kinematical characterization of the reaction products allows 
one to identify the particle emitting sources using invariant velocity plots 
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Figure 1: (Color online) Invariant velocity plot corresponding to all reaction products 
obtained in ^°Ca+^^C at 25 MeV/nucleon. The considered events have > 3. The 
vertical arrows indicate the projectile rapidity and its half. 
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Figure 2: (Color online) Yields of correlated a — a (a) and a-d (b) emissions out of QP 
expressed as a function of total kinetic energy (a) or excitation energy (b) and correspond- 
ing correlation functions (c and d). The insets in (a) and (c) correspond to zooms on the 
®Be(g.s.) peak. The inset in (d) details the a-d correlation function in the domain of 
^Li(2.19 MeV). Peak fits using Breit-Wigner distributions are illustrated with solid lines, 
while centroids are pointed by arrows. Full and open symbols in the insets correspond 
to different ways of calculating the angle under which the detected particles were emitted 
(see text). 

in the ■j/3± vs. 1/2 log [(^E + p\\) / (^E — p\\)] plane. E, py, p± stand for the 
energy and momentum components along and, respectively, perpendicular to 
the beam axis, /3 = f/c is the reduced velocity and 7 = — The 

invariant velocity plot corresponding to events having an a-multiplicity of at 
least 3 is depicted in Fig. [1] It shows the dominant binary character of the 
collisions with the formation of a quasi-projectile (QP) and a quasi-target 
(QT) not studied in this experiment. The QP source may be easily isolated 
by requiring a rapidity of reaction products larger than ?/pj.oj/2, where Uproj 
is the projectile rapidity. Hereafter, we shall focus exclusively on QP decay 
products with = 3. 

With the aim of investigating the reliability of the energy calibration, 
we show in Fig. |2] the a — a (top panel-left) and a-d (bottom panel-left) 
correlated spectra as a function of total kinetic energy in the center-of-mass of 
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the two particles (CM), Etot = Xli Ef^^\ and, respectively, excitation energy 
{Ef,x = Etot — Q)- For a — a correlation we considered the QP event sample 
described above. For a-d, we have additionally set the deuteron multiplicity 
to at least 1. The solid lines correspond to the peak fits using a Breit-Wigner 
distribution and peak centroids are pointed out by arrows. One may see that 
the a — a spectrum shows a narrow peak centered at 92 keV (r=84 keV) and 
a much broader peak centered around 3 MeV. The first peak corresponds to 
the ground state of *Be {Q=-92 keV) with r'^^^=5.57 eV and the second one 
to the first excited state at 3.03 MeV (r'^^^=1.5 MeV). The a-d spectrum 
shows a well-formed peak centered at 2.23 MeV (r=570 keV) corresponding 
to the first excited state of ^Li at 2.186 MeV (r^^P=24 keV). In both cases, 
the low energy levels are accurately determined (within 40-50 keV), meaning 
that the charged particle energy calibration is trustworthy. At variance with 
this, information on higher excited levels is blurred by particles originating 
from other decays whose abundance increases with energy. Nevertheless, 
contamination of spectra may be removed to a large extent by employing 
correlation function techniques which account for how much the correlation 
within the physical event differs from the underlying single-particle phase 
space. 

Correlation functions (CF) are defined as the ratio between the corre- 
lated (physical) yield Ycorr and the product of single particle yields, gener- 
ically termed as uncorrected spectrum Yuncorr, measured under the same 
conditions, 

1 + RiX) = j'-^f^y (1) 
J- uncorr\^ ) 

Yuncorr cau also bc built by mixing particles from different events as it is 
done in this work. In multi-particle CF 17|], the generic variable X is repre- 
sented by the total kinetic energy of the particles of interest in their center- 
of-mass frame Etot, by the excitation energy of their emitting source/state, 
Eex = Etot — Q or, in the case of two-particle correlations, by the relative mo- 
mentum. In nuclear physics, CF have been exploited to access spectroscopic 



properties [18|, give space-time information taking advantage of proximity ef- 



fects induced by Coulomb repulsion, and emphasize any production of events 



or sub-events with specific rare partitions [17|, 119|, 120(. The flat shape of the 



a-d CF around 1+R(X)=1 at E^^ >3 MeV (see (d) of Fig. [2]) confirms that 
the uncorrelated spectra are under control. 

Both correlated spectra and CF manifest a peak broadening. This is a 
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genuine consequence of detector finite granularity and energy resolution. Nu- 
merical simulations confirm that experimental distributions are compatible 
with an average energy resolution of detected alpha particles Re of 2%. The 
finite granularity is also responsible for a certain imprecision in the determi- 
nation of the energy. To illustrate this, the insets in Fig. |2] (a), (c) and (d) 
confront the correlated spectra and CF obtained under different hypotheses 
on the angle under which the detected particles were emitted. In the first 
case (solid symbols) we considered that all the particles which hit a certain 
module have been emitted under the angle corresponding to the geometri- 
cal center of that module. An alternative solution is to attribute to each 
particle a random angle in the domain allowed by the geometrical extension 
of the detector (open symbols). For all considered cases, the two scenarios 
lead to similar results, though the first hypothesis is preferable as producing 
smoother distributions and no extra broadening. In what follows only results 
obtained under the first hypothesis will be presented. 

Information on the a-particle unstable excited states of ^^C populated 
by the "^^Ca+^^C at 25 MeV/nucleon reaction may be extracted from the 
Sa-CF. The CF plotted in Fig. [3] corresponds to the whole set of events 
with rua = 3 and an uncorrelated spectrum built by full event mixing {i. 
e. each particle belongs to a different physical event). The error bars are 
calculated considering only the statistical errors of the correlated spectrum. 
Errors on the uncorrelated spectra have been reduced to negligible values 
by increasing the number of uncorrelated events as compared to correlated 
events. The CF shows two peaks centered at £'ea;=7.61 MeV (r=0.33 MeV) 
and £'ex=9-64 MeV (r=1.14 MeV). The first peak corresponds to the Hoyle 
state (i?|f^=7.654 MeV, T^^p=S.5 eV), while the second one is due to the 
complex region of excitations, characterized by the strong £^1^^=9.64 MeV 
(r'=^P=34 keV), 3" state and by the broad £^"^^=1 0.3 MeV, 0+ state submerg- 



ing a possible 2"^ state at 9.7 MeV j2l|, |22|. As before, the peak centroids 



reproduce well the experimental excitation energy values, while a significant 
broadening comes from the non-ideal detector properties (granularity and 
energy resolution). The inset of Fig. [3] illustrates the correlated spectrum 
in the energy domain of the Hoyle state and allows to estimate the rela- 
tive amount of genuine 3a decays versus background events. In the domain 
7.37-7.97 MeV this ratio amounts to 1. 

We demonstrated so far that the ^°Ca+^^C nuclear reaction at 25 MeV/nucleon 
populates excited states of ^^C nuclei which decay by 3-a emission. We now 
search for possible direct decays of the Hoyle state into equal energy parti- 
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Figure 3: (Color online) 3-a correlation function as a function of excitation energy. The 
arrows correspond to centroids of Breit-Wigner distributions (solid lines). Inset: zoom of 
the correlated spectrum in the energy domain of the Hoyle state. 



cles (DDE) and, eventually, estimate their branching ratio. The existence of 
DDE is manifest when looking at a two-dimensional CP in the (i?Q,)-RMS 
plane (Fig. IH^b)); (Ea) is the average kinetic energy of a-particles in their 
CM reference frame and RMS=-\^/ (£'^ ) — (En)'^. The uncorrelated spectrum 



is built by partial event mixing [23| of experimental events by taking two 
particles from the same event, in order to mimic for decay through ^Be, while 
the third one stems from a different event. The peak of the correlation func- 
tion localized at {E^) = 110 keV and very low RMS^ 25 keV corresponds to 
an equal sharing of the available energy of the Hoyle state among the three 
a particles. Note that 110 keV corresponds to the value of the maximum of 
the CP of Pig. 131 110 keV ~(7610+Q)/3. To go further in the interpretation 
of the CP and estimate branching ratios, the experimental results were com- 
pared to numerical Monte-Carlo simulations filtered by the multi-detector 
replica. These are performed by considering ^^C nuclei excited in the Hoyle 
state and boosted with a velocity distribution identical to the experimental 
one. They decay by one of the following mechanisms: (1) direct emission in 
three a-particles with equal energies (DDE), (2) sequential decay proceeding 
via the g.s. of '^Be with isotropic emission of the two a-particles from ^Be 
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(SD) and (3) direct a decay from a linear chain (DDL) [2J, |25| with an a 
at rest and an equal sharing of the available energy of the Hoyle state be- 
tween the two other a particles. Limitation of direct decays (DD), which 
in principle assume random energy sharing among the emitted particles, to 
the particular sub-class of DDE is due to the predicted condensate nature 
of the Hoyle state. In each case and, for sequential decay, at each decay 
step the available energy [E^.^ -|- Q) is shared among the reaction products 
so as to conserve the linear and angular momenta and, eventually, obey the 
postulated particularities of the decay. No preferential orientation between 
emitted particle velocities and source boost exists. The obtained reaction 
products are assumed to freely propagate toward the detector. Once the de- 
tector is reached, the events are filtered: if one or more particles hit a dead 
detector area or a detector which is out of work, the event is suppressed. 
If all particles have been detected, the velocity vectors are altered such as 
to account for finite angular and energy resolution. The simulated events 
are then analyzed as the experimental ones. The obtained detection effi- 
ciencies vary from 41 to 49% depending on the decay type. A comparison 
of filtered and normalized DDE (full), SD (dashed) and DDL (dot-dashed) 
simulations is depicted with thin lines in Fig. H] as a function of (a) 
and RMS (c). Very little sensitivity on the decay mechanism of Ycorr{,^E(^^ 
is observed, which is due to the angular resolution of the detection modules. 
By contrast, the kinetic energy dispersion (RMS) in the emitter CM man- 
ifests measurable sensitivity to the decay mechanism. For the CHIMERA 
granularity and a perfect energy resolution y'corr(-RAf 5) are peaked at 10, 70 
and 90 keV for DDE, SD and, respectively, DDL while for an average resolu- 
tion Re=2% the corresponding values are 30, 70 and 90 keV. This suggests 
that searching for the best agreement between experimental and simulated 
Ycorr{{Ea) , RM S) coustitutcs a pertinent procedure to quantify each decay 
channel. Not surprisingly, the general category of DD with random energy 
sharing (not considered hereafter) leads to a broad Ycorr{RM S) distribution 
which, under our experimental conditions, looks very similar to the domi- 
nating SD [l2|. For a realistic reproduction of the experimental data, we 
add to the pure simulated events the same proportion of background events 



that exists in the data. These are produced by partial event mixing [23 
of experimental events. The comparison is restricted to the energy domain, 
7.4 < < 7.8 MeV, to reduce the importance of background events to 
40%. Applying a minimization procedure and correcting for efficiencies, 
we infer that 7.5±4.0% of events correspond to DDE, 9.5±4.0% to DDL and 
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Figure 4: (Color online) Experimental three-a correlation function (b) expressed as a 
function of average kinetic energy - RMS of a particles corresponding to 1072 experimental 
events with 7.37 < E^x < 7.97 MeV. The uncorrelated yield is built such as to allow for 
decay through *Be. The dotted line marks the maximum RMS compatible with momentum 
conservation. (Ea) (a) and RMS (c) spectra (normalized to 1072 events) of simulated DDE 
(solid lines), SD (dashed lines) and DDL (dot-dashed lines) decays of the Hoyle state after 
filtering through the detector replica with Re— 2%. Panel (c) presents also the RMS 
projection of Ycorri{Ea) , RM S) /Yuncorr{{Ea) , RM S) (thick lines). 
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Figure 5: (Color online) RMS spectra of experimental (solid points) and best simu- 
lated events (solid line) in the region 90 < (Ea) < 150 keV. Experimental error bars are 
statistical. Dashed and dotted histograms illustrate simulation results if only sequential 
decay (60%) and, respectively, sequential decay (55.5%) + DDE (4.5%) are considered. In 
all cases, background (BG) events amount to 40%. 
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83.0±5.0% to SD. Error bars are estimated by taking into account statistical, 
and background errors. The RMS(i?Q) spectrum of simulated events (solid 
line) corresponding to the best and the energy domain 90 < (-Eq) < 150 
keV is displayed in Fig. |5] and compared to the corresponding experimental 
data (solid points). Note that a limitation to DDE and SD as possible decays 
(dotted histogram) leads to a significant increase of xLin without changing 
the DDE percentage. Results of simulations in which one allows only for 
SD are plotted with dashed lines for the sake of completeness. As one may 
notice, they correspond to a worse agreement with the data, particularly a 
sizeable underestimation of the number of low- RMS events. In all cases, the 
amount of background events is 40%. 

Finally we can now fully understand the CF pattern (Fig. Sl^b)); panel (c) 
of the figure presents also the simulated behaviours of Ycorr (RMS) / Yuncorr (RMS) 
(thick lines). Apart from DDE, the broad region around (Ea) = 90-130 keV 
and centered at RMS~ 70 keV corresponds to the sharing of the available 
energy between the two a's of ^Be and the remaining a of 191 keV. The 
peak at 130 keV and RMS = 90keV corresponds to the direct decay 

of a linear chain. Though the CF shows, as expected, a certain sensitivity to 
the event selection or event mixing recipe, its overall pattern remains stable. 
Particularly the DDL peak of the CF is systematically shifted by about 20 
keV with respect to the DDE one. As early mentioned, a small shift of the 
correlated spectrum is attributed to the finite size of detection modules (see 
Fig. m^a)). An extra shift could be introduced by statistical effects. 

A more popular way to visualize competing 3-particle decay mechanisms 
is the Dalitz plot. In this representation, DDE events must concentrate 



around the origin. Figure 1 in reference [26[ shows that our data manifest 
this pattern. Note also that, because of limited statistics and background, 
the rest of the figure cannot be interpreted. 

The same analyses have been performed for the complex region centered 
at 9.64 MeV where the statistics is much higher; no indication in favour of 
a direct 3-a decay with equal energies has been obtained |26|. For the Og 



state at 15.097 MeV of ^^O the statistics is too poor to allow unambiguous 
interpretation of low RMS events. 

In conclusion, the nuclear reaction ^°Ca-|-^^C at 25 MeV/nucleon bom- 
barding energy was used to produce states theoretically predicted as a- 
particle condensate states. Supposing that equal values of kinetic energy 
of the emitted a-particles represent a sufficient criterion for establishing the 
existence of a-particle condensation, we found that 7.5±4.0% of events cor- 
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responding to the Hoyle state decay fulfil this criterion. To our knowledge, 
this is the first direct experimental indication of a-particle condensation in 
nuclei. This study also evidenced the presence of direct alpha decays from 



a linear a-chain as recently theoretically mentioned [2J, |25[. An experiment 
with higher statistics is planned to study the ^^O case. 
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